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Results  of  this  measurement  program  on  and  around  Dulles  Airport  show: 

. Concorde  CO  emissions  from  taxiing  aircraft  dilute  to  background  levels  within 
2,000  ft  and  do  not  reach  the  terminal  in  measurable  amounts. 

. Emissions  from  the  airport  property  could  not  be  detected  at  Sterling  Park,  or 
several  more  distant  communities,  even  when  the  winds  were  blowing  toward  them 
from  the  airport. 

. Actual  Concorde  operations  ere  less  polluting  than  had  been  indicated  in  the 
Final  Environmental  Impact  Statement  (FEIS). 
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APPENDIX  A 


A. 1 Overview 

This  appendix  describes  in  detail  the  air  quality  and  supporting 
meteorological  data  obtained  at  Dulles  International  Airport  (1AD) 
during  intermittant  periods  between  15  May  1976  and  27  July  1977.  These 
data  were  acquired  as  part  of  the  monitoring  program  established  by 
the  Federal  Aviation  Administration  in  response  to  the  requirements  for 
air  quality  and  noise  monitoring  during  the  period  of  test  flights  by 
the  Concorde  into  Dulles  International  Airport. 

The  air  quality  data  consists  of  the  following  data  sets: 

• Hourly  avefaged  background  (May,  1976-September,  1976) 

• "Single  Event"  aircraft  CO  emissions  during  taxi  (May,  1976- 
August,  1976) 

• "Single  Event"  aircraft  NO^  emissions  during  takeoff  (September, 
1 976- April , 1977) 

• Multipoint  single-tower  measurements  of  CO  emissions  for  taxi 
(November,  1976) 

• Multipoint  two-tower  measurements  of  CO  emissions  for  taxi 
(February,  1976-April,  1977) 

• Multipoint  three-tower  measurements  of  CO  emissions  for  taxi 
(June-July  1977) 

Supporting  meteorological  data  (other  than  wind  speed  and  wind  direction 
obtained  on-site,  temperature  and  temperature  gradient  for  the  two-tower 
and  three-tower  experiment)  include  the  following  sets: 

• Hourly  surface  observations  as  reported  by  NOAA  at  the  airport 

• Twice  daily  rawinsonde  ascents  made  at  the  NOAA  upper  air 
station  near  the  airport. 

Additional  supporting  air  quality  data  were  obtained  from  the  network 
of  monitoring  stations  established  by  agencies  of  the  states  of  Virginia 
and  Maryland,  and  by  the  National  Institute  of  Health.  These  data  are 
included  in  the  hourly-averaged  background  data  set. 
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All  of  the  data  were  obtained  either  as  strip  chart  recordings  or 
as  written  log  entries.  The  strip  charts  were  subsequently  digitized, 
and  the  information  formatted,  coded  and  put  on  computer  punch  cards 
to  facilitate  machine  analysis  and  processing.  In  the  following  sections 
of  this  document,  the  content  and  format  of  these  data  sets  are  discussed. 
The  locations  of  the  monitoring  sites  in  the  vicinity  of  the  airport, 
and  the  characteristics  of  these  sites  are  presented  in  Figure  A-l  and 
Table  A-l. 

A .2  Hourly  Averaged  Background  Data  Set 
A. 2. 1 Station  Set 

This  set  of  data  consists  primarily  of  measurements  from 

e Site  *1  - Sterling  Park 
• Site  *6  - South  Rump 

and  the  stations  from  the  states  of  Virginia  and  Maryland,  and  the 
National  Institute  of  Health.  These  stations  are  identified  as  follows: 


State 


Virginia 


Maryland 


Station  Name 

Herndon 

Massey 

Xerox  Training  School 
Club  Run 
Lewinsville 

Douglas  Elementary  School 
Seven  Corners 

Burning  Tree 

Gaithersburg  High  School 
Poolesville  Elementary  School 
Rockville 

National  Institute  of  Health  ■ 
Bethesda 


The  location  of  these  regional  stations,  relative  to  Dulles  International 
Airport  is  shown  in  Figure  A- 2.  Table  A-2  summarizes  the  parameters 
measured  at  each  site.  The  units  used  to  specify  the  parameters  are 
identified  in  Table  A-3. 
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Figure  A-l  Air  Monitoring  Sites  - Dulles  Airport 


TABLI-  A-l 


CHARACTERISTICS  OF  MONITORING  SITES 


Slta  # 

Approximate 

Location 

Maaaurauant 

Function 

Operational 

Statue 

1 

Starling  Park 

CO,  NO  4 NO., 
total  HC,  par- 
ticles, osone 
aathana,  WD/WS* 

Monitor  air  quality  change 

In  the  coaaunlty  during 
Concorde  Operation 

June-Sept  (1974) 

2 

2000  ft  north 

4 123  ft.  aaat 
of  tha  aaat 
runway** 

CO,  NO*.  WD/WS* 

Measure  tekaoff  salsslone  ae 

precursor  to  defining  loca- 
tions for  3-statlon  takeoff 
grouping,  background  a tat ion 

May-Aug  (1974) 

3 

1100  ft.  north 

4 1000  ft  aaat 
of  tha  a~  it 
runway 

MO, 

Measure  takeoff  emissions 
as  precursor  to  defining 
locations  for  3-statlon 
takeoff  grouping. 

May-Aug  (1974) 

4 

Taxi  Grouping  - 
480  f b north 
of  tha  Jet 
ranp  taxiway  4 
200  ft. aaat  of 
tha  aoblla 
lounga  rasp 

CO,  WD/WS* 

Trace  emission  propagation 
during  taxi  (slngla  event) 

May-Sept  (1974) 

S 

Taxi  Grouping  - 
190  ft. north  of 
tha  Jat  raap 
taxiway  4 . 200 
ft. aaat  of  tha 
■oh lla  lounga 
raap 

CO 

Seas  aa  4 

May-Sept  (1974) 

4 

South  adga  of 
tha  aaln  raap; 
on  tha  aouth 

accaaa  road 

00,  NO  4 NO 2, 
total  HC,  par- 
tlclaa,  oaona 
WD/WS* 

Monitoring  pollution  back- 
ground of  tha  airport 

May-Sept  (1974) 

7 

Northwoat  of 
airport  (M44- 
ownod) 

Vertical  teap. 

Moaeure  Inver a Ion  baas 

Continuous 

• 

Want  of  tha 
want  runway 
(WOAA-ownad) 

WD/WS* 

(no  racordar) 

Monitor  wind  apaad  and 
direction 

Continuous 

» 

3000  ft.  aouth 
of  Slta  4 

(MQAA-ownod) 

WD/WS* 

direction10*  *****  *** 

Continuous 

10 

Taxi  Grouping  - 
Midway  batwaan 
Sltaa  443 

00 

Sane  aa  4 

June-July  (1974) 

*WD/US-Ulnd  direction /Wind  speed 


**A11  dimensions  ara  measured  fro*  cantarllna  of  rasp,  runway  or  taxivav 
unlaaa  otharvlaa  naadad. 
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^-1  (Cont  inued) 


Site  # 

Approximate 

Location 

Measurement 

Function 

Operational 

Status 

XI 

Taxi  Grouping  - 
200  f t. north  of 
Slta  4 

CO 

Same  ae  4 

July-Sept  (1976) 

12 

Takeoff  Group- 
ing - 285  ft 
east  and  100 
ft. north  of 

Site  13 

CO,  NO  6 N02, 
total  HC,  ozone, 
WD/US* 

Trace  emission  propagation 
during  takeoff  (single  event) 

Oct  76- May  1976 

13 

Takeoff  Group- 
1 lng  - 185  ft. 
east  & 140  ft. 
north  of  Site 

14 

CO,  NOx 

Same  as  12 

Sept  76- May  1976 

14 

Takeoff  Grouping 
450  ft. east  of 
the  east  runway 

6 1040  ft. south 
of  Its  north  end 

CO,  NOx 

Same  as  12 

Sept  76-May  1976 

15 

16 

17 

Start/ldle 
Grouping  - north 
of  the  west  & of 
the  Jet  reap 
taxiway 

CO,  WD/WS*, 

(at  one  site) 

Trace  emission  propagation 
during  engine  start/ldle 
(single  event) 

Spot  check 

18 

South  of  the 
east  runway 

NOx 

Monitor  takeoff  emissions 

Spot  check 

19 

North  of  the 
west  runway 

NO* 

Monitor  landing  emissions 

Spot  check 

20T 

South  edge  of 
main  ramp  1700 
ftt**  west  of 
Runway  19L, 

56-ft. elevation 
on  tower 

Air  Intake  position  (tower) 

November  (1976) 

2 IT 

41-ft.  elevation 
on  tower 

- 

Sams  as  20T 

November  (1976) 

22T 

26-ft.  elevation 
on  tower 

Same  ae  20T 

November  (1976) 

TABLE  A-l  (Continued) 


Site  I ' 

Approximate 

Location 

Measurement 

Function 

Operational 

Status 

2 ST 

14-f t. elevation 
on  tower 

- 

Same  as  20T 

November  (1976) 

24 

South  edge  of 
main  ramp  ** 

1665  ft*** 
west  of  Runway 

19L 

CO,  VS/VD* 

Tower  measurements 

November  (1976) 

25 

1665  ft***  weat 
of  Runway  19L 

164  f t south 
kf  Site  24 

CO 

Tower  measurements 

November  (1976) 

26 

11700  ft***  west 
of  Runway  19L 

164  ft. south  of 
Site  25 

Air  intake  position 
(Surface) 

November  (1976) 

27T 

Same  as  20T 

Soft  ele- 
vation on 
tower 

Same  as  20T 

Feb-March  (1977) 

26TT 

164  ft. south 
of  20T.  ao-ft. 
elevation  on 
tower 

Same  as  20T 

Feb-Aprll  (1977) 

29  TT 

Same  as  28TT 

56  -ft.  ele- 
vation on 

tower 

Same  as  20T 

Feb- April  (1977) 

SOTT 

41-ft. 

elevation  on 
tower 

• 

Same  as  20T 

j 

Feb-Aprll  (1977) 

3 ITT 

26-f  t.  ele- 
vation on 
tower 

Same  as  20V 

Feb- April  (1977) 

32TT 

14 -ft.  ele- 
vation on 

1 

Same  as  20T 

Ftb- April  (1977) 

i 

tover 


***215  ft.  south  of  south  )et  ramp  centerline 
****South  end  centerline 
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Figure  A- 2 Regional  Air  Quality  Monitoring  Stations 


BACKGROUND  STATIONS 


*TSP  values  are  not  reported  on  an  hourly  basis  but  are  reported  as  daily  (24)  hour  values. 


TABLE  A- 3 


UNIT  SPECIFICATIONS  FOR  BACKGROUND  STATION  DATA  SET 

Parameter  Units 


Sterling  Park  § 

South  Ramp 

Virginia 

Maryland 

Wind  Direction  (WD) 

Degrees  (N=360°) 

Wind  Speed  (WS) 

MPH 

Carbon  Monoxide  (CO) 

ppm 

mg/m 

7 

ppm 

Nitrogen  Dioxidp  (NO^) 

ppm 

yg/ni 

7 

ppm 

Nitric  Oxide  (NO) 

ppm 

yg/m 

*7 

ppm 

Ozone  (03) 

ppm 

yg/m 

none 

Total  Hydrocarbons  (THC) 

ppm 

mg/m 

7 

none 

7 

Total  Suspended  Particulates  pg/m3 

Pg/m3 

yg/ni 

yg/jn 

A. 2. 2 Data  Formats 


A. 2. 2.1  Sterling  Park  and  South  Ramp 

The  outputs  of  the  sensors  were  recorded  by  means  of  strip  chart 
recorders  running  at  a speed  of  1 inch  per  hour.  These  charts  were 
returned  to  ERT' s Concord  facilities  where  they  were  digitized  auto- 
matically. The  output  of  the  digitizer  was  fed  into  a computer  where  a 
tape  of  hourly- averaged  values,  in  the  proper  units,  was  created  for 
each  of  the  parameters.  The  card- image  format  of  the  final  data  set  as 
shown  in  Table  A-4.  Examples  of  the  information  content  of  the  back- 
ground data  set  are  shown  in  Figures  A-3  and  A-4.  In  Figure  A-3,  the 
data  have  been  stratified  by  site,  parameter  (NO.,),  month  (August  and 
September),  day  of  month,  and  hour  of  day.  In  Figure  A-3,  the  strati- 
fication is  by  stability  class  and  wind  direction.  With  the  parameters 
contained  in  the  data  set,  other  kinds  of  stratified  analyses  are 
possible. 


A. 2. 2. 2 State  Agency  Data  Set 

These  data,  received  by  ERT  through  the  FAA,  consists  of  hand-entry 
tabulation  sheets  in  the  standard  SAROAD  format  used  by  the  Environmental 
Protection  Agency  (EPA) . An  example  of  this  entry  form  is  shown  in 
Figure  A-S.  The  format  is  already  set  up  for  computer  processing.  The 
data  were  put  on  computer  punch  cards  following  the  SAROAD  format  using 
the  parameter  specifications  as  entered  in  the  data  sheets.  A data 
tape  was  prepared  for  this  data  set. 

A. 2. 2. 3 National  Institute  of  Health 

The  data  received  from  the  National  Institute  of  Health  consisted 
of  computer  printouts  of  hourly-uvcraged  values  similar  in  format  to 
that  shown  as  Figure  A-3.  An  example  of  the  print  out  is  shown  in 
Figure  A-(>.  The  computer  tabulated  data  were  put  on  computer  punch 
cards  in  a format  identical  to  that  used  for  the  Sterling  Park  and 
South  Ramp  data  set. 
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TABLE  A-4 


EXTENDED  ENVIROMAP  ARCHIVE  FORMAT 


Field  Cols  Variable  iFormat  I Description 


I FORM 


IliXP  (1) 

(2) 

(3) 

(4) 

(5) 

(6) 


I YEAR 
1 DAY 


I HOUR 
1 MINUTE 


20  INCNT 

21 -24  KODE(l) 

2S-34  VAL(l) 

3S-37  KCT(l) 

38  KFG(l) 

39- -10 


Dataset  Indcntification 


Experiment  ID 
Experiment  Mode 
Data  Type 
Interval  Units 
Integrat  i on  • I ntorval 
Process  O|)tion 


Year  (YY) 
Day  (ODD) 


Hour  (1111) 
Minute  (MM) 


Data  Item  Count 


A4  Sensor  ID  Code:  1st  Sensor 

I PE  10. 3 Data  Value  . 

0P13  Sample  Count 
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Figure  A- 4 txample  of  Background  Data  Stratified  by  Stability 


Figure  A-5  SAAOAD  Format 
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Figure  A-6  Example  of  the  NIH  Data 


A. 3 Single  Event  and  Multi-Point  Data  Set  Description 


A. 3.1  Measuring  and  Recording  Techniques 

The  data  from  these  series  of  experiments  consist  of  the  following: 

a.  High  speed  (2  cm/min)  strip  chart  recordings  of  the  outputs 
of  each  sensor 

b.  Observer  log  of  events  and  the  instantaneous  wind  speed  and 
direction  at  the  time  of  an  event,  where  an  event  is  defined 
as  the  taxi  of  an  aircraft  in  front  of  the  line  of  CO  monitors 

Figure  A-7  shows  an  example  of  the  typical  3-pen  trace  recorded  on  the 
high  speed  chart  recorder  for  an  event. 

The  following  information  was  digitized  from  the  strip  charts  for 
each  sensor  output  and  for  each  event: 

• Background  Concentration  (ppm) 

• Peak  Concentration  (above  background)  (ppm) 

• Event  pulse  half  width  (secs) 

• Peak  time  (secs) 

• 1/2  peak  time  (secs) 

The  definition  of  these  parameters  are  shown  schematically  in  Figure  A-8 
Note  that  the  half  width  time  and  peak  time  is  measured  relative  to  the 
event  time.  Event  time  is  the  time  recorded  by  the  observer  in  the 
observer  log  as  the  time  that  the  aircraft  visually  passes  in  front  of 
the  CO  monitoring  shelters.  The  recording  of  these  times  are  not  of 
sufficient  precision  to  be  used  to  determine  the  transport  time  of  the 
plane  from  the  source  to  the  monitoring.  The  times  proivde  a measure 
of  reconstructing  the  pulse  shape  should  this  be  a desirable  parameter. 
The  above  set  of  parameters  together  with  observer  records  of 

• Date 

• Event  time  (LST) 

• Aircraft  type 


Concentration 


O 


“ 

Figure  A-8  Pollution  Parameter  Definition 
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• Direction  of  aircraft  movement 

• Aircraft  operation  mode 

• Wind  speed  and  direction 

form  the  basis  of  the  data  card  of  an  event. 

A general  computer  card  format  was  designed  as  shown  in  Figure  A-9 
to  incorporate  the  digitized  strip  chart  values  of  CO,  N0x  and  THC, 
and  the  field  observer  logged  information  pertaining  to  the  passage 
of  each  aircraft.  The  following  is  a description  of  the  data  identifica- 
tion fields  used  on  the  punched  cards: 


Card  Column 

1.  2 


3-10 


11 


12 


13 


Parameter 

Code  or 

Name 

Range 

Description  and  Comment 

Event  No. 

01-99, 

Event  100  is  coded  as  00, 

00-nn . 

however  they  are  numbered 
chronologically  within  each 
experiment  series. 

Time  Period 

MO 

Month 

DY 

Lay 

HR 

Hour  (1-24) 

MN 

Min 

Aircraft  Type 

1 

B-707 

2 

B-727 

3 

B-747 

4 

DC-8 

5 

L1011 

6 

Concorde 

7 

Turbo 

8 

DC- 10 

8 

B737  (three  tower  only) 

9 

DC  9 

0 (zero) 

DC-10  (three  tower  only) 

b (blank) 

Cessna  (three  tower  only) 

Mode  of  Opera- 

1 

ENGINE  start 

tion 

2 

Take  off 

3 

Landing 

4 

Taxi 

Direction  of 

1 

East  ) 

Operation 

2 

West  1 

Omitted  in 

3 

North  1 

one-tower  set 

4 

South  J 
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Card  Column 

Parameter 

Name 

Code  or 
Range 

Description  and  Comment 

14,  15 

Instantaneous 
Wind  Speed 

00-98 

Miles  per  hour 

Note:  99  is  missing  value 
indicator 

16,  17,  18 

Instantaneous 
Wind  Direction 

50.0-99.0 

Percent  scale  on  strip 
chart  recorder  (100%  = 540° 
50%  = 0°) 

Note:  99.9  is  missing  value 
indicator 

19,  20 

50  meter  timing 

0.0-9. 8 

(tenths  of  seconds) 

(omitted  in  one-tower  set) 
Note:  9.9  is  missing  value 
indicator 


The  first  20  card  columns  are  repeated  within  each  set  of  cards  corres- 
ponding to  a specific  event  for  each  experiment  series.  Any  exceptions 
appear  in  the  above  'Description  and  Comments'  column.  In  figure  A-9, 
the  computer  card  format  is  shown.  The  fields  from  card  columns  21 
through  80  contain  the  digitized  strip  chart  values  of  air  quality  data, 
and  are  described  by  sensor  location  within  the  next  two  subsections  of 
Singe  Event  and  Multipoint  experiment  descriptions. 

A. 3. 2 Single  Event  Experiment  Description 

Outputs  of  3 CO  sensors  arranged  in  a N-S  direction  to  the  north  of 
the  taxiway  to  Runway  19L  and  3 NO^  sensors  arranged  in  a SW-NE  line  to 
the  east  side  of  the  end  of  Runway  19L  constitute  the  data  base  for  the 
two  "Single  Event"  experiments.  Since  only  3 sensors  were  operative 
during  each  experiment,  the  digitized  values  obtained  from  the  strip 
chart  recorders  are  contained  on  a single  card  per  event. 


A. 3. 3 Tower  Demonstration  Data 

A. 3. 3.1  Experiment  Description 

Three  series  of  experiments  were  conducted  using  meteorological 
towers.  The  first  made  use  of  a single  tower  located  near  the  south 
edge  of  the  main  ramp,  west  of  Runway  191..  Two  additional  "ground 
level"  CO  sensors  were  arranged  inline  southward  from  the  tower,  involving 
a total  of  6 CO  sensors.  The  ground  level  tower  sensor  also  recorded 


■X 
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THC.  The  second  experiment  consisted  of  the  addition  of  1 CO  sensor  at 
the  top  of  the  first  tower,  and  the  replacement  of  the  middle  ground 
level  sensor  with  a second  identical  tower.  With  the  third  ground  level 
sensor,  a total  of  11  CO  sensors  were  involved.  This  two  tower  experi- 
ment grouping  also  included  the  measurement  of  THC  at  the  ground  level 
position  of  the  first  tower  as  well  as  the  addition  at  the  same  location 
of  wind  speed,  wind  direction,  temperature  and  temperature  differential. 

For  the  final  three  tower  series  of  experiments  the  most  distant 
ground  level  sensor  was  replaced  by  another  tower.  This  provided  an 
additional  two  CO  sensors,  located  at  the  top  and  middle  positions 
corresponding  to  the  same  first  and  second  tower  sensor  locations.  A 
total  of  13  CO  sensors  were  involved.  An  additional  wind  speed  and  wind 
direction  sensor  was  placed  at  the  top  of  the  first  tower.  Lower  level 
meteorological  data  and  THC  continued  to  be  monitored  at  the  same  loca- 
tions as  used  during  the  Pouble  Tower  experiment,  with  the  exception 
that  the  lowest  sensor  levels  were  at  six  feet  instead  of  fourteen  feet. 
Important  note:  during  the  Double  Tower  series,  wind  speed  and  wind 
direction  data  were  recorded  on  the  same  chart  with  scale  of  0-50®  = 
wind  speed,  and  50  to  100®  full  chart  scale  representing  wind  direction. 
During  the  Triple  Tower  experiment  series,  the  recorder  pens  were 
reversed,  and  wind  direction  is  recorded  from  0-50%  and  wind  speed  from 
50  to  100%.  [Both  scale  values  are  the  same.  One  chart  unit  WS  = 

1 mph,  1 chart  unit  WD  = 10.8°,  full  scale  WS  = 50  mph  full  scale 
WD  = 540°.] 

A. 3. 3. 2 Data  Set  Description 

Three  separate  data  sets  were  generated  for  each  experiment  group- 
ing. This  is  due  to  the  unique  combination  of  parameters  and  their 
increase  in  quantity  of  sensors  as  the  tower  experiments  progressed. 

There  are,  however,  common  elements  to  each  tower  series  data  set  which 
are  incorporated  into  the  computer  card  format  used  in  the  basis  of  the 
identification  of  an  event.  These  pertain  to  the  field  observers 
records  of  aircraft  passage,  as  described  above.  The  sensor  fields  are 
valid  for  the  Double  and  Triple  Tower  series,  however  the  Single  Tower 
data  omitted  the  peak  and  half-peak  tine  values.  The  following  describes 
the  air  quality  and  meteorological  data  as  coded  on  computer  punch 
cards : 
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1)  ONE  TOWER  DATA  - CO  (one  card/event) 


Card  Column 

Sensor  ID 

Parameter 

21-23 

24-26 

20T 

-divisions  above  background  (F3.1) 
-pulse  width  @ half  peak  (sec)  (F3.0) 

27-29 

30-32 

21T 

Repeated  as  for  20T 

33-35 

36-38 

22T 

Repeated  as  for  20T 

39-41 

42-44 

23T 

Repeated  as  for  20T 

45-47 

48-50 

25G 

Repeated  as  for  20T 

51-53 

54-56 

26G 

Repeated  as  for  20T 

57 

none 

Blank 

58-60 

Background 

2)  TWO  TOWER 

DATA  (four 

cards  per  event) 

a.  CO  data  - 

Since,  as  shown  in  Figure  A-9,  there  are  four  allowable 
sensors  per  card,  the  following  comprise  the  data  set: 

Card  1:  Sensors  27T,  20T,  21T,  22T. 

Card  2:  Sensors  23T,  28TT,  29TT,  30TT. 

Card  3:  Sensors  31TT,  32TT,  26G,  23THC. 

b.  Met  data  - 

Card  4 - The  first  20  card  columns  are  identical  to  the 
first  three  cards.  The  remaining  columns  contain  measured 
and  digitized  wind  speed,  direction  and  turbulence  data 
according  to  the  following  format: 
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Card  Columns 
21-25 

26-30 

31-35 

36-40 

41-45 

46-50 

51-55 

56-60 


Parameters 

3 min.  ave.  wind  speed,  u,  (mph) , for 
time  centered  upon  event 

3 min.  ave.  wind  direction,  0,  (deg) 
for  time  centered  upon  event 

°g  (deg)  for  3 min.  sample  of  6-sec 
averagea  values 

(mph)  = Og  • u • (1  rad/57.3  deg) 

AT  = temperature  difference  and  first 
tower  between  14  and  67  ft.  3 min 
average  starting  at  clock  interval 
closest  to  event  time 

AT^  = same  as  AT,  but  18°  min  average 
centered  on  events  3 min.  period 

T ambient  air  temperature,  hourly 
average  estimate. 

Pasquill  Turner  stability  class  (4  = neutral) 


c.  THREE  TOWER  DATA  SET  - CO  data  (4  cards  per  event) 

Card  1:  Sensors  27T,  20T,  21T,  22T 
Card  2:  Sensors  23T,  28TT,  29TT,  30TT 
Card  3:  Sensors  31TT,  32TT,  34TTT,  33TTT 
Card  4:  Sensors  26TTT,  23THC 


A. 3. 3. 3 Strip-Chart  Data  Summary  for  Tower  Experiments 

Available  are  4 volumes  containing  xerox  reductions  of  strip  chart 
images  categorized  by  aircraft  event,  separated  by  date  corresponding  vo 
each  series  of  tower  experiments.  Table  A-5  is  a summary  by  aircraft 
types.  Tables  A-6(a)  - A-6(c)  provide  a distribution  of  events  by  date 
for  the  Single,  Double  and  Triple  Tower  experiment  series.  Contained 
within  each  volume  is  a listing  of  each  event  number,  time  and  type  of 
aircraft  passage. 
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TABLE  A- 5 


SUMMARY  BY  AIRCRAFT  TYPE 
l-OR 


ONE,  TWO 

AND  THREE  TOWER 

EXPERIMENTS 

Aircraft  Type 

One  lower 

Two  Tower 

Three  Tower 

Tot  a 1 

B 707 

4t> 

31 

49 

126 

B 727 

15 

34 

27 

7o 

B 747 

15 

3 

8 

20 

0C- 8 

11 

5 

12 

28 

L 1011 

17 

0 

15 

32 

SST 

9 

1 

9 

19 

Turbo 

0 

0 

1 

1 

B 737 

0 

0 

5 

5 

nc-r 

2 

0 

3 

s 

DC- 10 

21 

0 

'18 

39 

Cessna 

0 

0 

3 

3 

A -37 

0 

0 

:i 

21 

Tot  a 1 

ISO 

74 

173 

383 

Table  A-6  a 

DISTRIBUTION  OP  EVENTS  BY  DATE 
I Single  Tower 


Date 

Event  No. 

11/1/76 

1 - 

15 

11/4/76 

16  - 

20 

11/5/76 

21  - 

33 

11/7/76 

34  - 

45 

11/8/76 

46  - 

59 

11/10/76 

60  - 

86 

11/11/76 

87  - 

99 

11/12/76 

100  - 

123 

11/13/76 

124  - 

137 

11/15/76 

(2  unnumbered) 
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Table  A-6  b 

DISTRIBUTION  OP  EVENTS  BY  DATE 


l 


II  Double  Tower 


Date  Event  No. 


1 

2/21/77 

1 - 4 

2 

2/26/77 

5 (6  f,  7) 

3 

2/28/77 

8-12 

4 

3/2/77 

13  - 19 

5 

3/4/77 

20  - 32 

6 

3/7/77 

33  - 45 

7 

3/8/77 

46  - 50 

8 

3/15/77 

51  - 62 

9 

3/19/77 

63  - 66 

10 

3/23/77 

67 

11 

3/24/77 

68,69,70 

12 

3/2b/77 

71  - 78 

13 

3/31/77 

79,80 
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Table  A-6  c 


: 


DISTRIBUTION  OF  EVENTS  BY  DATE 
III  Triple  Tower 


Date  Event  No. 


1 

6/27/77 

1 - 18 

2 

7/6/77 

19  - 23 

3 

7/14/77 

24  - 50 

4 

7/15/77 

51  - 58 

5 

7/22/77 

59  - 75 

6 

7/26/77 

76  - 129 

7 

7/27/77 

130  - 100 

8 

7/72/77 

161  - 174  (A- 37  only) 
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B.  SINGLE  EVENT  IMPACT  MODELING 


B.l  Statistical/Empirical  Models 

The  advantage  of  statistical  models  is  that  they  summarize  the 
observational  data  in  an  informative  manner  as  well  as  suggesting  the 
relative  importance  of  measured  variables,  without  relying  upon  elaborate 
knowledge  of  the  mechanisms  producing  the  results.  Statistical  methods 
are  most  efficiently  employed  when  a specific  hypothesis  is  first 
formulated  for  testing.  When  statistical  methods  are  used  to  "explore" 
the  data  for  meaningful  relations,  there  is  always  a finite  chance  that 
a "statistically  significant"  correlation  or  regression  will  be  found 
where  none  actually  exists.  It  is  thus  best  to  examine  the  physical 
relationships  of  the  variable  chosen  for  correlation  or  regression 
analysis  so  that  highly  cross-correlated  variables  can  be  excluded,  to 
avoid  confusing  the  analysis. 

The  single  event  twin-tower  experiments  have  been  summarized, 
(Appendix  C)  stratified  according  to  aircraft  type,  to  compare  the 
means,  standard  deviations,  and  ranges  of  the  most  important  variables. 
Simple  linear  correlations  and  cross  plotting  of  measured  concentrations 
and  estimated  doses  vs.  wind  and  turbulence  variables  and  aircraft  speed 
are  carried  out  to  provide  the  initial  insights  into  the  strengths  of 

I 

expected  trends  or  relationships.  The  multilinear  regression  process 
makes  efficient  use  of  each  measuremcn;  by  allowing  it  to  contribute 
information  about  the  effects  of  several  variables  simultaneously.  As 
in  Section  a relation  of  the  form: 


E(y)  * *0  + B1  X1  * 62  x2  * •••  * 0m  xm 

where  E(y)  is  the  expected  value  of  the  dependent  variable  y,  and  Xj , 

x,...x  are  the  values  of  the  m independent  variables,  is  sought. 
l m 

For  the  development  of  a single  event  model,  the  dependence  of  con- 
centrations and  doses  upon  the  wind  speed  and  direction  and  turbulence 
intensities  examined  initially.  As  expected,  the  plume  transport  path 
geometry  and  factors  affecting  plume  rise  were  important  when  a sufficient 
number  of  cases  were  examined.  Since  aircraft  type  is  a variable  that 
does  not  have  a continuous  scale  of  values,  stratification  of  the  data 
set  by  that  variable  helped  to  maintain  regression  model  sensitivity. 
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B.2  Transport  Models 


As  discussed  in  Section  7,  pollutant  transport  and  dispersion  models 
which  attempt  to  describe  the  pathway  and/or  concentration  distribution 
of  materials  downwind  of  their  source  are  frequently  employed  as 
alternatives  to  statistical  interpretations  of  experimental  results. 

Both  the  Gaussian  models  and  numerical  advection  diffusion  calculations, 
described  below,  derive  from  the  differential  equation  describing  the 
conservation  of  mass  under  advection  and  diffusion: 


3_ 

3t 


♦ u 


3x 


♦ v 


3y 


♦ w 


3z 


x = k v x 


(B-2) 


where  x is  the  pollutant  concentration  and  u,  v and  w are  the  components 
of  the  mean  wind  velocity  in  the  x,  y and  z directions,  respectively. 

The  Gaussian  form  of  the  solution  was  seen  by  Gifford  (1960)  to  be 
analogous  to  the  empirical  diffusion  equations  used  by  Sutton  (1932)  and 
Pasquill  (1959),  and  it  is  that  form  that  has  come  into  most  frequent 
use.  As  presented  in  the  Turner's  Workbook  of  Atmospheric  Dispersion 
Estimates  (1970),  the  general  form  of  the  equation  for  estimating  pollu- 
tant concentrations  from  a continuous  point  source  is: 


x(x,y,z,H) 


2no  o u 
y z 


g—  exp  (-1/2  (-i)2] 


[exp  [-1/2  * exp  [-1/2  t5-— )2]1 

w w 7 J 


(B-3) 


where 


(x,y,z)  are  the  (along  wind,  cross-wind,  and  vertical)  components 
of  a Cartesian  coordinate  system 

X is  the  pollutant  concentration  (mass/ volume) 

H is  the  effective  height  (stack  height  plus  plume  rise)  of 

emission,  and  therefore  the  centerline  height  of  the 
plume  (length) 

Q is  the  source  strength  (mass/time) 
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y' 


a 
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are  dispersion  coefficients  that  are  measures  of  cross- 
wind  and  vertical  plume  spread.  These  two  parameters  are 
functions  of  downwind  distance  and  atmospheric  stability 
(length) 


u is  average  wind  speed  (length/time) 


The  source  base  is  at  z = 0 in  the  coordinate  system,  and  the  plume 
centerline  reaches  the  equilibrium  height  H (plume  rise  plus  stack 
height)  at  some  distance  downwind  from  the  source.  The  most  important 
assumptions  upon  which  the  equation  is  based  are  the  following: 

1)  The  wind  speed  and  direction  are  constant  throughout  the 
period  of  interest. 

2)  The  plume  rises  until  it  reaches  an  equilibrium  altitude; 
thereafter  the  plume  centerline  height  remains  constant  at  all 
further  downwind  distances. 

3)  The  distribution  of  concentration  values  off  the  centerline  is 
given  by  the  product  of  two  Gaussian  distributions,  ore  in  the 
y-direction  and  one  in  the  z-direction. 

4)  The  concentration  profiles  described  by  the  Gaussian  form  are 
not  "instantaneous"  plume  profiles;  they  represent  concentra- 
tions averaged  over  one  hour.  Consequently,  they  incorporate 
the  normal  variability  of  wind  flow  for  this  time  period. 

5)  None  of  the  effluent  is  lost  from  the  plume.  Therefore,  when 
the  plume  intersects  the  ground  surface,  it  is  assumed  that 
all  material  is  reflected  back  above  the  ground. 

6)  The  effluent  rate  is  constant,  and  the  meteorological  param- 
eters determining  plume  geometry  are  constant;  (i.e.,  the 
equation  represents  steady  state  conditions). 

Because  this  form  of  the  Gaussian  model  is  not  entirely  appropriate 
for  emissions  from  a single  aircraft  traveling  down  a taxiway  or  a 
runway,  the  next  three  sections  describe  each  of  three  modified,  or 
quasi-instantaneous,  approaches  to  this  Gaussian  form  of  model.  Also 
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described  is  a numerical  advection  method  for  solving  the  same  original 
transport  equation.  All  of  these  approaches  are  expected  to  be  more 
appropriate  than  equation  (B-3)  for  modeling  the  concentrations  in  the 
near  field  of  a taxiway  or  runway. 

B.2.1  Quasi- Instantaneous  Gaussian  Models 

When  single-event  concentrations  are  to  be  predicted,  the  instan- 
taneous line  source  form  appears  immediately  to  be  the  most  correct 
form,  since  we  are  concerned  with  the  passage  of  a single  source  (the 
airplane)  along  a line  (the  runway).  When  many  airplanes  are  in  con- 
tinuous transit  along  the  runway  with  only  negligible  spacing  in  between, 
it  is  normal  to  sum  up  their  contributions  and  arrive  at  the  usual 
"continuous  line  source"  form  of  the  Gaussian  equation  commonly  used  in 
highway  modeling.  In  the  analysis  of  the  single-event  data  (Section  6), 
three  different  modeling  methods  have  been  compared  against  measurement 
data  obtained  at  Dulles  Airport. 

These  methods  differ  only  in  the  derivation  of  plume  spread  para- 
meters and  o,.  In  each  case,  the  derived  parameters  are  utilized 
with  a single  "instantaneous  line  source"  form  of  the  Gaussian  diffusion 
equation  to  determine  the  expected  peak  concentration  at  each  sensor. 

The  instantaneous  line  source  equation  is  derived  below. 

Instantaneous  Line  Source  Gaussian  Equation 


In  order  to  derive  the  Gaussian  diffusion  equation  for  an  instan- 
taneous line  source  at  height,  H,  it  is  convenient  to  rearrange 
Equation  B-3: 


x(x,y.z,H) 


2 _L 

u /Tiro 


. I X_ 


1 ilM 

2 o 


(B-4) 


where  x is  the  pollutant  concentration  in  g/m  , Q is  the  source  strength 
in  g/sec,  u is  the  mean  wind  speed  in  m/sec,  and  and  o,  are  the 
standard  deviations  of  plume  read  in  the  lateral  and  vertical  direc- 
tions, respectively. 


This  form  of  Gaussian  aquation  allows  for  "perfect  reflection"  of 
pollutants  by  folding  over  the  vertical  distribution  at  ground  level. 

The  result  is  a doubling  of  the  surface  concentration  through  the  second 
of  the  curly  bracket  terms,  as  can  be  shown  by  setting  the  height,  z,  to 
zero. 

The  two  curly  bracketed  terms  are,  respectively,  the  "probability 
densities"  of  pollutant  mass  in  the  y and  z directions.  The  area  under 
each  probability  density  curve  is,  by  definition,  unity. 

The  conversion  from  the  continuous  point  source  form  to  the  instan- 
taneous line  source  form  of  the  Gaussian  equation  is  straightforward. 

The  continuous  point  source  distribution  is  Gaussian  along  the  y-axis, 
and  constant  along  the  x-axis.  Conversely,  for  the  instantaneous  line 
source,  the  distribution  is  Gaussian  along  the  x-axis  (normal-to-the- 
runway  axis),  and  constant  along  the  y-axis  (runway  axis).  The  proba- 
bility density  along  the  x-axis  may  be  written  as 


P(x) 


(B-S) 


where  Xm  is  the  distance  along  the  normal  from  the  plume  centerline  or 
parallel  to  the  mean  wind  component  normal  to  the  runway.  Since  the 
distance,  Xm,  is  measured  along  the  normal  wind  axis  from  the  plume 
centerline,  we  have  to  take  into  account  the  motion  of  the  plume  center- 
line  as  advected  by  the  wind  normal  to  the  runway.  If  x is  measured  in 
a stationary  coordinate  system  from  the  runway  centerline  to  a sensor  at 
distance  d1(  then  the  transformation  to  the  moving  coordinate  system 
(Xm)  yields: 


u t 
n 


(B-b) 


1 


Where  un  is  again  the  component  of  the  mean  wind  normal  to  the  runway. 

When  0 is  the  angle  between  the  mean  wind  direction  and  the  normal  to 

the  source  line,  then  u = ucos  0.  The  coordinate  system  outlined  above 

n 

is  depicted  in  Figure  (B-l).  The  probability  density  along  the  x-axis 
now  coincides  with  that  given  by  Sutton  (1932),  as  used  by  Turner  (1970) 
in  equation  S.21: 
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Here,  P(x)  reaches  a maximum  when  d,  = unt,  or  the  centerline  of  the 
plume  blows  over  a sensor  at  distance,  d7,  from  the  runway  centerline. 

The  only  remaining  problem  is  that  of  defining  the  proper  line 
source  emission  term,  g [gm/m] . The  term  Q/u  in  the  continuous  point 
source  equation  does  have  the  proper  units,  but  we  must  redefine  the 
velocity  u in  the  transformation  to  the  instantaneous  line  source  form. 
Three  velocity  vector  components  along  the  runway  axis  can  serve  to 
stretch  and  compress  the  plume  and  alter  the  emission  rate  per  unit 
meter.  They  are:  (1)  the  component  of  the  mean  wind  speed,  u^  = u ‘sine, 
parallel  to  the  runway  (along  the  y-axis);  (2)  the  jet  exhaust  velocity 
Uj  (y) ; and  (3)  the  aircraft  taxi  speed,  s. 

Lacking  satisfactory  knowledge  of  the  jet  exhaust  velocity  as  a 
function  of  down-runway  distance,  only  the  effects  of  s and  Up  have  been 
considered  in  this  analysis;  so  that 


q 

rather  than 


-..3— 

s + u 

p 


q = Q/(s  + u.  (y)  ♦ Up) 


(B-8) 


(B-9) 


Substituting  P(x)  (Equation  B-7)  for 

2 

exp  (-1/2  ^~)  = P(y)  (B-10) 

/JrT  o y 

y 

in  Equation  B-3  and  using  the  redefined  source  term  (B-8),  the  resultant 
form  of  the  instantaneous  line  source  equation,  as  used  in  the  preliminary 
analysis  is: 
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In  trial  analyses,  plume  rise  was  characterized  by  three  different 
methods.  The  first  characterization  assumes  no  plume  rise,  so  that  the 
mean  plume  height  stays  at  the  mean  exhaust  emission  height.  The  other 
two  characterizations  are  the  expressions  fit  statistically  by  R.  J. 
Yammartino  (March  1S77)  to  the  data  from  the  November  1976  single-tower 
experiment.  His  first  expression  is  based  on  a mean  plume  height,  the 
second  on  a Gaussian  plume  height. 

At  the  first  tower,  the  three  plume  rise  characterizations  express 
the  plume  centerline  height  as 


H = Al,2  + Bl,2//J7 


where  u(mph)  is  the  wind  speed  for  each  single  event.  The  subscripts  1 
and  2 refer  to  two  sets  of  aircraft  types  (see  Table  7-2).  Values  for 
the  parameters  are: 


No  Plume  Rise  Mean  Height  Gaussian  Centerline 


Aj(ft)  mean  exhaust  height  24.7 
A* (ft)  mean  exhaust  height  29.7 
Bj  (ft  mph  ''i')  0 20.-5 
B*(ft  mph*'  j 0 9.9 


-2.6 

9.1 

65 

43 


In  the  analysis,  plume  rise  was  assumed  to  be  complete  by  the  time 
the  plume  reaches  the  first  tower.  Thus,  H was  assumed  to  be  identical 
at  all  sensor  distances  downwind  of  the  taxiway  or  runway. 

Modeling  Approaches  Used  to  Define  Sigmas 

In  order  to  evaluate  the  Gaussian  diffusion  equation,  the  standard 
deviations  of  plume  growth,  and  a ^ , must  be  known.  This  analysis 
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explores  the  methods  used  to  quantify  and  o through  three  methods: 

(1)  a method  suggested  by  Turner  using  Pasquil 1-Gifford-Tumer  disper- 
sion parameters,  (2)  a method  based  on  the  measured  turbulence  data  at 
Dulles  airport  during  the  experiment,  and  (3)  a circular  jet  method. 

The  first  two  methods  are  chosen  because  of  their  frequent  usefulness 
in  obtaining  estimates  of  concentrations  or  doses  from  Gaussian  models. 

The  3rd  method,  the  circular  jet  approach,  is  chosen  in  order  to 
realistically  model  the  physical  process  of  jet  engine  wake  dispersion 
(see  Lin,  1958) . 

At  time,  t = 0,  all  three  methods  assume  an  instantaneous  line 
source  centered  along  the  runway,  as  indicated  in  Figure  (B-l).  At  the 
aircraft  itself,  the  plume  is  assumed  to  have  an  initial  volume  based  on 
the  engine  geometry  of  the  particular  aircraft  being  modelled.  Behind 
the  aircraft,  the  volume  of  the  plume  grows  in  direct  proportion  to  the 
growth  rates  of  and  o^.  When  the  initial  size  of  the  plume  is  greater 
than  zero,  the  usual  course  is  to  describe  a virtual  source  distance, 

Sy,  at  which  the  plume  dimensions  extrapolate  to  a zero  source  volume. 

Thus,  ox  and  will  be  always  proportional  to  the  distance  from  the 
virtual  source. 

Dispersion  parameters  at  any  down-plume  distance  can  then  be  expressed 
as  functions: 


(B- 12) 


where  d^  is  the  distance  of  the  sensor  normal  to  the  runway  at  which 
ox  and  oz  are  to  be  quantified. 

From  the  plume's  initial  location  along  the  runway  at  time,  t = 0, 
the  plume  is  advected  by  the  component  of  the  mean  wind,  u^,  normal  to 
the  runway.  Figure  (B-l)  illustrates  this  process  as  the  plume  is 
advected  by  un  * ucosS  from  the  runway  to  the  sensor  at  dz  = 165.5m  = 
unt.  During  this  same  time  t,  the  portion  of  the  plume  at  the  sensor 
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has  travelled  a distance  from  the  plane,  parallel  to  the  runway,  d 


l 


(s  ♦ u^") t , where  u^  ■ usinO  is  the  component  of  the  mean  wind  parallel 
to  the  runway  (neglecting  the  contribution  of  t ho  jet  exhaust  velocity, 
u.(y)  in  transporting  plume  mass  away  From  the  aircraft).  Therefore, 
the  distance  that  exhaust  is  transported  from  the  aircraft  to  the 
sensor  is: 


|/,2  ,2  , I /.  s ♦ usini) 

I/d  ♦ d - d 1/1  ♦ — 

r 1 f ucosO 


(h- 13) 


so  that,  in  equation  (B-12)  above 
s - f (0  ♦ S ) 

X v V 

s,  - g (D  ♦ Sv)  (B-14) 


where  the  functions  f and  g depend  on  the  method  used.  A description  of 
each  method  and  the  critical  assumptions  necessary  in  their  use  follows. 

Pasqui 1 1 -Turner  Method 

The  initial  volume  of  the  plume  is  taken  so  that  4.72  o equals 
the  exhaust  span  and  4.72  equals  twice  the  height  of  the  highest 
engine  above  ground.  Thus,  the  initial  plume  width  in  each  dimension 
is  taken  as  two  pulse  widths  (one  pulse  width  at  half  peak  height  equals 
2.5oo).  This  is  the  same  as  saying  that  99%  of  the  plume  mass  is  con- 
tained between  the  rwo  most  widely  spaced  exhaust  jets  in  the  horizontal 
and  between  the  ground  surface  and  twice  the  height  of  the  highest 
exhaust  jet  in  the  vertical.  These  initial  source  dimensions  are 
described  pietorally  in  Figure  (B-l).  Observations  of  aircraft  exhaust 
plumes  closely  support  these  choices  of  initial  source  dimensions.  Two 
virtual  source  distances  are  described,  one  based  on  the  initial  at 
the  aircraft,  and  one  based  on  the  initial  o at  the  aircraft.  (See 
equation  B-12.)  The  sigmas  are  then  quantified  at  sensor  distances 
d,  ■ t>5,  115,  105  meters  by  the  look-up  tables  in  Turner  (1970,  Figures  3.2 
and  3.3). 
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Critical  assumptions  to  this  method  are: 


a)  The  growth  rate  of  the  standard  deviation  normal  to  the 
runway,  o*,  can  be  described  by  growth  curves  from  Turner 
(1970) . 

b)  Although  a quas i-instantancous  line  source  is  being  modeled, 
the  o ' s and  o 's  chosen  are  those  used  to  describe  spread 

X 2 

from  cont inuous  sources.  This  choice  is  based  upon  the  fact 
that  the  source  is  somewhat  more  persistent  than  a puff- like 
source . 

c)  The  only  contribution  of  the  jet  exhaust  dynamics  is  to  "advect" 
pollutant  mass  parallel  to  the  runway.  This  approach  assumes 
that  jet  exhaust  dynamics  plays  no  role  in  diffusing  the  plume 
normal  to  the  runway. 

d)  The  effect  of  wake  caused  by  the  body  of  the  aircraft,  as  it 
obstructs  the  moan  wind  flow,  can  be  entirely  taken  into 
consideration  by  the  choice  of  the  above  mentioned  initial 
source  volume. 

e)  At  the  runway  centerline,  the  plume  is  assumed  to  have  a 
Gaussian  distribution,  although  one  might  expect  the  plume  to 
be  composed  of  the  superimposed  Gaussian  distributions  of 
spread  from  each  individual  engine  exhaust. 

f)  The  plume  does  not  grow  appreciably  in  the  time  its  width  from 
♦2o  to  -2o  passes  the  sensor.  This  assumption  has  no  effect 
on  the  centerline  plume  prediction  but  will  cause  errors  in 
off-centerline  predictions  and  derived  pulse  lengths.  An 
allowance  for  plume  spread  during  transit  past  a sensor  should 
be  made  in  order  to  compare  predicted  pulse  lengths  to  obser- 
ved pulse  lengths.  Experimental  data  suggests  that  the  plume 
grows  appreciably  as  it  passes  a sensor;  many  of  the  pulses 

on  the  strip  chart  are  log-normal. 
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g)  The  taxi  path  of  the  aircraft  is  sufficiently  long  to  allow 
a steady-state  situation  to  develop.  The  model  will  be  in 
error  if  the  aircraft  turns  at  the  end  of  the  taxiway  while 
it  is  still  contributing  pollutant  mass  to  any  sensor.  Such 
cases  should  be  analyzed  in  a different  manner,  as  by  the 
Argonne  Line  Source  Model. 

Turbulence  Measurement  Data  Method 

The  standard  deviations  of  plume  growth  and  o^,  can  be  derived 
at  any  downwind  distance  by  knowing  either  the  standard  deviation  of 
the  azimuthal  wind  angle,  oQ,  or  the  turbulence  intensity,  cru/u. 

Based  on  o0,  the  relationships  are: 

o = ao  (D  + S 
x 6 v v * 
x 

°Z  = coe  (D  + sv  )d  (b-15) 

z 

where  a,  b,  c,  and  d are  coefficients  tabulated  in  Slade,  1965 
(page  134,  Table  4.7). 

Based  on  c>u/u,  the  relationships  are  (Pasquill,  1974): 

°x  = °u(D  + Sv  )/u  (B-16) 

x 

a = 0.6  a for  neutral  conditions 

z x 

= for  moderately  unstable  conditions 

Both  of  these  approaches  were  tested  in  the  analysis.  The  param- 
eters o.  and  o /u  were  derived  from  the  measured  turbulence  data,  and 
D+Sy  was  computed  as  in  the  Pasquill-Tumer  method.  Critical 

assujftp^fo^s  a)  and  c)  through  g)  from  the  Pasquill-Tumer  method  also 
apply  here. 
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Axisymmetric  Jet  Method 


It  is  well  documented  (see  Lin,  1958J  that  an  axisymmetric  (cir- 
cular) jet  in  a boundary-free  medium  growing  in  width,  De>  with  respect 
to  downwind  distance  can  be  described  by: 

2o  * D ■ 2d  tan  (B-17) 

x e 

where  d is  the  downplume  distance  from  the  source,  and  $ = 5°  is  the 
radial  growth  angle. 

The  modeling  method  here  is  still  quite  similar  to  the  Pasquill- 
Tumer  method  and  can  also  be  described  by  Figure  (B-l).  The  only  real 
difference  here  is  that  the  above  linear  growth  formula  is  used  to 
derive  instead  of  Turner's  sigma  curves.  Also,  since  the  growth  is 
circular  in  the  circular  jet  approach,  is  assumed. 

The  virtual  source  distance  Sv>  can  be  defined  through  trigonometry 
as: 


S 

v 


Hi 

2 tan  $ 


(B- 18) 


where  is  the  exhaust  span  in  the  horizontal. 

Only  one  virtual  source  distance  is  described  here  and  it  is 
dependent  only  on  the  horizontal  exhaust  span  and  independent  of  the 
height  of  the  exhaust  jets.  This  is  critical  to  the  circular  jet  method 
since  the  method  is  based  on  the  assumption  that  the  plume  always  has 
equal  horizontal  and  vertical  dimensions  and  it  is  not,  therefore, 
theoretically  correct  to  describe  the  plume  as  anything  but  circular. 
Thus,  the  horizontal  exhaust  span  depicted  for  the  Pasqui 11 -Turner 
and  Turbulence  Data  Measurement  Methods  in  Figure  (B-2)  is,  for  the 
Circular  Jet  Method,  also  taken  as  the  initial  vertical  width. 

The  complete  formula  for  and  o„  is: 


o 


x 


o 


z 


(B- 19) 
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The  concentration,  x»  is  computed  as  before.  The  critical  assumptions 
to  this  approach  are: 

a)  Diffusion  is  a consequence  of  only  the  dynamics  of  the  jet 
exhaust  plume  itself.  The  contribution  of  mixing  due  to 
turbulence  in  the  ambient  air  has  been  considered  small,  in 
this  preliminary  analysis.  In  later  analysis,  ambient 
turbulence  may  be  shown  to  become  important  after  some  down- 
plume  distance  where  jet-generated  turbulence  has  decayed  to 
such  an  extent  as  to  have  comparable  intensity  or  eddy 
diffusity  as  the  situation  in  the  ambient  air. 

b)  A basic  assumption  is  that  we  are,  indeed,  describing  a 
circular  jet,  even  though  the  horizontal  exhaust  span,  L^,  is 
much  longer  than  the  vertical  exhaust  span,  Lv<  Indeed,  the 
approximation  of  a plane  jet  (very  long  width,  very  short 
height)  may  be  more  valid.  In  this  case,  Lin  (1958)  suggests 
using  a half-width  growth  angle  in  the  horizontal  of  6-1/2° 
rather  than  5°.  The  half-width  growth  angle  in  the  vertical, 
however,  is  not  well  documented,  since  few  experiments  can 
sufficiently  match  the  plane  jet  in  its  theoretical  limit: 
infinite  width,  zero  height. 

Assumptions  (e) , (f)  and  (g)  from  the  Pasquil  1-Tumer  method 
above  also  apply  to  this  approach. 

B.2.2  Quasi-Continuous  Gaussian  Models 

The  transport  model  utilized  in  the  AVAP  program  is  based  upon  the 
Argonne  Puff  Line  Model  approximation  to  the  continuous  Gaussian  form  of 
line  source  model  identified  at  the  beginning  of  the  previous  section. 
The  Puff  Line  Model  is  conventionally  used  in  AVAP  with  source  and 
meteorological  data  and  dispersion  parameters  that  most  nearly  represent 
hourly  average  behavior  of  a continuous  source.  In  the  current  model 
development  effort,  the  basic  Puff  Line  Model  has  been  modified  to 
incorporate  "instantaneous"  (three  minutes)  meteorological  data  and 
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dispersion  parameters.  The  assumption  of  a continuous  rate  of  emission 
over  the  length  of  the  source  aircraft's  path  is  made  here,  but  the 
assumption  is  also  made  that  the  period  of  emission  is  short  compared  to 
the  transport  time  to  the  sensor  receptors.  A more  approximate  assump- 
tion that  accompanies  use  of  this  form  for  analysis  is  that  the  period 
of  peak  concentration  passage  (or  dose  accumulation)  is  small  compared 
with  the  (three  minute)  averaging  time  for  the  meteorology.  The  advan- 
tage of  the  Gaussian  Puff  Line  Model  is  its  ability  to  accurately 
calculate  dispersion  rates  for  wind  angles  which  are  nearly  parallel  to 
the  taxiway,  and  the  ability  to  handle  accelerating  or  decelerating 
sources.  A detailed  description  follows.* 

To  treat  long  thin  sources  and,  in  particular,  mobile  sources  such 
as  aircraft  and  motor  vehicles,  a general  line  source  model  has  been 
developed.  The  basic  theory  of  the  Line  Source  Model  is  presented  in 
this  section.  For  further  discussion  of  the  theory  the  interested 
reader  is  referred  to  (Wang  and  Rote,  1975;  Rote  and  Wangen,  1975). 

The  basic  line  source  equations  can  be  derived  using  an  approach 
similar  to  that  used  in  the  Puff  Model.  We  begin  by  assuming  that  the 
effluent  emitted  over  a time  duration  i,  from  a finite  straight  line 
segment,  can  be  treated  as  a sequence  of  long  thin  "puffs"  or  "linear 
puffs"  extending  over  the  length  of  the  line  segment.  In  other  words, 
if  the  average  concentration  at  the  receptor  is  given  by: 


X 


C(t)dt, 


(B-20) 


where  C(t)  is  the  time  dependent  concentration,  then  we  are  assuming 
that  this  expression  can  be  replaced  by  the  alternative  expression: 


X x l \i  = [ 7 | C.(t)dt  (B-21) 

i i J 

where  C^(t)  represents  the  contribution  to  the  receptor  concentration 
from  the  ith  linear  puff  as  a function  of  the  time.  The  duration  x is 
taken  to  be  the  averaging  time  (=  1 hour)  over  which  the  meteorological 
parameters  are  considered  constant.  It  is  further  assumed  that  the 

•This  description  has  been  prepared  by  the  original  authors  of  the  Puff 
Line  Model,  D.  M.  Rote  and  L.  E.  Wangen. 
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time  of  formation  of  each  puff  is  short  compared  to  the  averaging 
time  t and  that  each  puff  rapidly  comes  to  rest  relative  to  the  ambient 
air  mass.  The  effect  of  the  penetration  of  a jet  exhaust  plume  into  the 
ambient  air  due  to  the  high  exhaust  velocity  is  discussed  in  Appendix  B. 

With  these  assumptions,  it  is  convenient  to  treat  the  transport 
of  a linear  puff  of  pollutant  using  the  Green's  function  technique.  For 
a ground- level  horizontal  line  source  at  an  angle  4>  relative  to  the 
wind,  we  can  write  the  concentration  due  the  release  of  such  a puff, 
at  the  receptor  point  (x,y,z)  averaged  over  the  period  t as: 


,L 

d K K (x,y,z;£,t)  q^. 
'o 


(B-22) 


where  q^  is  the  linear  mass  density  of  the  puff,  and  L is  the  length 
of  the  straight  line  segment. 

In  general,  the  linear  mass  density  q will  be  a function  of  position 
along  a particular  line  source.  For  the  present,  only  uniform  linear 
mass  densities  will  be  considered.  (The  topic  of  non-uniform  emission 
densities  will  be  discussed  in  detail  later.)  For  the  special  case  of 
constant  density, 


r ' 

Mass 

i 

Length 

J 
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where  = the  total  mass  of  the  puff.  K (x,y,z;£,t)  is  known  as  the 
transport  kernel  and  can  be  expressed  as 


K (x,y,z;£,t)  = (2it)  " oL  o 


- . -3/2  -2  -1 
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exp 
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where  $ is  the  angle  of  the  line  relative  to  the  x-axis.  In  writing 
down  the  above  expression,  we  have  temporarily  ignored  ground  reflection 
and  lid  effects  and  chosen  the  x-axis  of  the  coordinate  system  along 
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the  wind  vector  y with  the  origin  of  the  coordinates  at  one  end  of  the 
line.  Furthermore,  we  have  assumed  that  the  turbulent  diffusion  of 
pollutant  about  the  centerline  drawn  downwind  from  each  infinitesimal 
line  element  d£  can  be  represented  in  terms  of  the  usual  Gaussian 
horizontal  and  vertical  dispersion  coefficients  and  o z>  respectively. 
These  coefficients  are  assumed  to  be  functions  of  the  stability  and 
travel  time. 

If  t = the  travel  time  for  a puff,  then,  since  we  have  assumed 
that  t ^ is  short  compared  to  t,  it  follows  that  the  exposure  time  of 
the  receptor  to  the  puff  will  be  limited  to  the  time  interval  from 
t - s^/2  to  t + s^/2,  where  the  exposure  time  duration  (s^)  is  less 
than  x.  In  other  words,  the  transport  kernel  K (x,y,z;C*t)  is  highly 
peaked  in  the  neighborhood  of  t so  that  one  can  effectively  rewrite 
the  integral  in  Equation  B-22  as 


dt 


d£  qi  K (x,y,z;£,t) 


(B-24) 


Furthermore,  provided  K (x  y,z;£,t)  is  sufficiently  peaked  at  t = t,  it 
is  permissible  to  replace  the  functions  (t)  and  a ^ (t)  by  their  values 
at  t * t,  i.e.,  (t) , a„  = a (t)  inside  the  integral.  Consequently, 

for  uniform  emission  density,  the  average  concentration  can  be  approximated 
by 


X,  = 2 (2ir) 
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The  division  of  a long  line  into  short  segments  to  insure  the 
criterion  s.  < t is  performed  automatically  within  the  line  source 
dispersion  model  computer  program.  However,  no  attempt  is  made  within 
the  computer  program  to  subdivide  the  width  of  the  line.  Hence,  for  a 
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line  of  width  W,  after  a sufficiently  long  travel  time  T,  the  condition 
> r will  be  violated.  In  other  words,  assuming  that  the  downwind 
plume  dimension  is  given  by  2 * (J,  T ♦ T^)  where  J is  the  stability, 

T ♦ Tw  is  the  total  travel  time  from  a pseudo  upwind  zero  width  line 

n 

source,  a value  of  T exists  such  that,  for  greater  travel  times. 


oh  (J,  T ♦ Tw)  > x * WS. 


(B- 26) 


For  the  extreme  conditions  of  stability  J * 1,  t « 1 hour,  and  WS  = 1 m/sec, 
we  find  that,  using  the  dispersion  curves  of  Figure  13,  s^  > t when 


T ♦ T > 3,000  sec.  Under  such  conditions  an  approximation  to  the 

n 


hourly  average  concentration  could  be  computed  for  any  given  hour 
during  the  exposure  of  the  receptor  to  the  plume  by  multiplying  by 

T/S.. 

When  the  line  is  elevated  to  height  H and  ground  reflection  is 
taken  into  consideration  the  z dependent  factor 


exp 
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should  be  replaced  by 
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When  one  end  of  the  line  is  at  the  origin  and  the  other  extends  essentially 
out  to  infinity,  the  first  error  function  in  Equation  (B-25)  goes  to  1, 
and  Equation  (B-25)  reduces  to  the  form 
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(B-28) 
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The  dispersion  coefficients  oy  (J,T)  and  o,  (J  ,T)  already  defined  in 
connection  with  the  short-term  point  and  area  source  model  are  also  used 
for  o.  and  o,,  respectively.  The  initial  width  Ay  and  height  At  of  a 
physical  line  source  are  treated  in  analogy  with  the  physical  point 
and  area  sources  by  first  assigning  initial  values  of  horizontal  and 
vertical  dispersion  o and  o,  and  then  computing  the  corresponding 
positions  of  pseudo  upwind  line  sources.  Returning  to  the  total 
contribution  due  to  a train  of  linear  puffs,  one  can  regard  the  train 
of  puffs  as  resulting  from  a sequence  of  discrete  events  such  as  the 
movement  of  a number  of  aircraft  or  in  the  case  of  a continuous  stream  of 
automotive  traffic  as  simply  a mathematical  artifice.  In  the  former 
case,  to  compute  the  total  x.  one  can  simply  multiply  by  the  number 
of  aircraft  movements  provided  that  they  occur  within  the  appropriate 
time  period.  In  the  latter  case,  Q.  is  replaced  by  the  total  emission 
Q for  the  line  segment  for  the  time  period. 

When  the  angle  between  the  wind  and  the  line  becomes  small  and 
the  line  is  short  (see  below),  liquation  B-28  can  be  approximated  by 
its  limiting  expression  at  = 0, 
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The  general  expression  for  an  inclined  line  source  with  angle  of  inclina- 
tion 0 is  derived  by  Rote  and  Wangen,  1975. 

In  applying  the  above  formulas  to  actual  situations,  it  is  important 
to  bear  in  mind  that  the  line  geometry  must  be  such  that  an  instantaneous 
release  from  each  line  segment  i must  contribute  to  the  receptor 
significantly  only  during  an  interval  small  compared  to  t,  i.e.,  s^  < t. 
Thus,  for  a long  line  source  at  small  angles  relative  to  the  wind  one 
should  divide  the  line  into  shorter  segments  and  evaluate  effective 

and  for  each  segment  separately.  It  is  precisely  these  considera- 
tions that  prevent  one  from  writing  down  a simple  expression  such  as 
Equation  B-25  for  the  long  line  in  the  small  angle  case.  If  a line 
segment  is  above  the  lid  or  if  an  inclined  line  penetrates  the  lid,  that 
portion  above  the  lid  is  excluded  by  the  calculation. 
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B.3  Dispersion  Parameter  Averaging  Times 


P 

» 


According  to  Pasquill  (1974)  the  relationship  between  the  appro- 
priate averaging  time  for  wind  and  turbulence  parameters  and  plume  puff 
travel  time  is: 

S » = */u6 

where  S is  the  averaging  time,  T is  the  puff  travel  time  to  the 
receptor  (monitor)  at  distance  X,  and  8 is  the  ratio  of  the  Eulerian 
(measured  at  a fixed  point)  to  the  Lagrangian  scale  (moving  with  the 
mean  wind  u).*  As  an  example,  for  a puff  moving  at  u*  3 m/sec  (^7  mph) 
the  appropriate  averaging  time  for  measuring  the  lateral  speed  of  a 
puff  at  a monitor  60  m away  is  (60/4.5)  = 5 sec  as  long  as  the  sampling 
(or  emissions)  time  is  at  least  five  times  the  20  sec  travel  time. 
Selection  of  the  optimum  sampling  and  averaging  times  for  standard 
reporting  of  wind  parameters  is  complicated  however,  by  four  factors; 

(1)  the  times  are  functions  of  wind  speed,  (2)  they  are  functions  of 
monitors  distance  from  the  source,  (3)  the  time  period  of  release  is 
difficult  to  determine  (it  may  range  from  "vl  to  30  sec  for  moving  air- 
craft, while  it  may  be  virtually  continuous  in  a cueing  situation), 

(4)  the  use  of  chart  speeds  of  1-2  cm/min  makes  it  difficult  to  accurately 
digitize  data  for  averaging  periods  less  than  6 sec,  and  finally, 

(5)  the  value  of  6 is  not  known  to  have  a constant  value  for  all  atmo- 
spheric stabilities.**  Consideration  of  all  of  these  factors,  and  the 
relatively  large  data  base  of  research  results  that  exist  for  sampling 
times  of  3 min  leads  to  the  conclusion  that  fixed  sampling  and  averaging 
times  should  be  used  for  the  estimation  of  turbulence  parameters  in  the 
present  analysis  program.  The  present  analysis  results  were  based  on  the 
use  of  high  speed  wind  charts  will  have  wind  direction  and  speed 
digitized  every  6 sec  for  the  3 min  period  centered  upon  the  passing  of 

a monitored  aircraft.  This  will  enable  calculation  of  o • x = o and 

o v 

/u  • x 3 o*.  These  values  can  be  used  to  estimate  o,.  the  vertical 
dispersion  parameter  from; 


•Pasquill  has  established  that  6=4  ,s  a good  estimate  by  its  mean 
value . 

**The  Hay  Pasquill  measurements  of  6 cover  a range  of  l.l  to  8.5. 
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These  estimates  of  lateral  dispersion  parameters  will  be  directly 
comparable  to  the  Pasqui 11 -Turner  estimate  for  the  50  to  200  m 
distances  over  which  continuous  concentrations  monitored.  For  greater 
distances  these  parameters  should  optimally  be  calculated  according  to 
the  most  recent  recommendations  of  Pasquill  (1976). 

B.4  MINUIT*  Program  Used  in  Multiparameter  Fitting  Models 

The  CERN  program  MINUIT  was  used  in  the  parameterization  model. 

This  program  will  minimize  a function  of  N variables,  compute  the 

2 2 

covariance  matrix,  and  find  at  x = + 1 true  errors.  Three  types 

of  minimization  methods  are  available  to  the  user.  They  may  be  briefly 
characterized  as  follows: 

1)  Search  for  the  minimum  using  a Monte  Carlo  method.  This  may 
be  used  at  the  beginning  of  a fit  when  no  reasonable  starting 
point  is  known,  or  when  it  is  suspected  that  there  are  several 
minima.  It  does  not  converge  in  the  usual  sense. 

2)  Minimize  using  the  simplex  method  of  Nelder  and  Mead.**  This 
is  a "safe"  and  reasonably  fast  method  when  far  from  the 
minimum.  It  does  not  compute  the  covariance  matrix,  but  it 
gives  some  estimate  of  its  diagonal  elements. 

3)  Minimize  using  a variable  metric  method  by  Davidon.t  This  is 
extremely  fast  near  a minimum  or  in  any  "nearly-quadratic"  region 
but  slower  if  the  function  is  badly  behaved.  It  fails  com- 
pletely in  regions  where  the  covariance  matrix  is  not  positive- 
definite. 


*F . James  and  M.  Roos,  MINUIT,  CERN  D506,  D516  1971. 

**J.  A.  Nelder  and  R.  Mead,  The  Computer  Journal,  7,  308  (1967). 
+W.  C.  Davidon,  The  Computer  Journal,  _10,  406  (1968). 
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STATISTICAL  SUFtlARY  OF  INFIELD  CO  TESTS  FOR  DC  10 
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Correlation  of  Concentrations  with  Meteorological 
Variables:  B 707  Aircraft/14  ft  Level  CO  Measurements 
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Correlation  of  Concentrations  with  Meteorological 
Variables:  B 707  Aircraft/41  ft  Level  CO  Measurements 
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Figure  D1.31  Diurnal  Analysis 
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Figure  D1.3S  Diurnal  Analysis 
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